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Node-Reordering Method for Stiffness Matrix 
Wavefront Reduction in Structural Analysis 


It is anticipated that in the performance of struc- 
tural analysis, the solution of structural matrices by 
means of the present bandwidth concepts will be re- 
placed by the more efficient wavefront concepts. A 
method has been developed which fills an important 
void by providing a readily-implemented approach 
to automatic node relabeling that is consistent with 
requirements of the wavefront concept. Specific ap- 
plications are in the analysis of aircraft, building 
structures, radar and surveillance structures, bridges, 


etc., or any other structure that is studied with the 


aid of a large and complex analytical model. 



The purpose of reordering for an analytical model 
of a particular structure is to minimize the number of 
degrees of freedom that determine the maximum of 
the bandwidth or the wavefront of the stiffness matrix. 
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Nodal Connectivity Matrix 

Frequently, the number of degrees of freedom per 
node is approximately the same for all nodes. There- 
fore, it is reasonable to perform reordering on the 
basis of the nodes and nodal connectivity, which is 
considerably simpler than on the basis of degrees of 
freedom and stiffness matrix connectivity. A new 
algorithm, described here, has been found effective 
in reducing the wavefront 
Wavefront counting for a sample structure is illus- 
trated in the diagram. The wavefront at a row of the 
nodal connectivity matrix is the number of active 
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columns that follow the diagonal element; a column 
becomes active at the row for which there is the first 
entry for that column and it remains active until that 
column is absorbed into the diagonal of the connec- 
tivity matrix. In the nodal connectivity matrix, heavy 
vertical lines are drawn for each column beginning 
at the row when the column becomes active and end- 
ing at the diagonal. To obtain the wavefront at each 
row, it is necessary only to count the number of these 
heavy lines crossing each row. 

Alternatively, the wavefront at any row can be 
computed as equal to the wavefront at the preceding 
row, plus the number of new columns that become 
active. If the column corresponding to the row is 
currently in the wavefront, then the wavefront is re- 
duced by one. The subtraction is made to account for 
columns that leave the wavefront by virtue of reach- 
ing the diagonal. The isolated column to the right of 
the connectivity matrix in the diagram gives the 
wavefront at the row. The maximum wavefront is 
equal to 6, which occurs at both the second and third 
rows; this can be reduced by application of the re- 
sequencing procedure. 

The recommended resequencing approach can be 
called a “minimum growth” method; that is, assuming 
that resequencing has already been used to identify 
the nodes that will constitute the new first n nodes, 
then the node selected as the new (n+l)th node is 
the one for which the wavefront at the (n+l)th node 
will represent the smallest increase with respect to 
the wavefront at the nth node. In accordance with the 
alternative method described for wavefront counting, 
the increase can be a positive integer, zero, or neg- 
ative unity. Very often, any one of several nodes will 
give the same minimum wavefront increase for the 
(n+l)th node; when this happens, the first node is 
selected, for this tends to make the new resequencing 
as close as possible to the original. Any one of several 
methods can be used to choose the node to be first 
in the new sequence. 

The operations for resequencing can be conven- 
iently organized in a tableau format that employs the 
original nodal connectivity matrix and a new table 
that is constructed, to contain the change in wavefront 
for each possible selection of the next node to be 
placed in the new sequence. Tableau resequencing 


can be performed without difficulty for matrices with 
orders of from 20 to 30; when the order becomes 
much larger, the tableau becomes cumbersome and 
it is advisable to perform resequencing by computer. 

A computer program has been written to automate 
the procedure by performing essentially the same 
operations as in executing the tableau. The largest 
model for which the program has been used to date 
has a connectivity matrix with over 500 nodes. In one 
test of the program for this model, the initial se- 
quencing had been developed by inspection and 
study of the connectivity. The maximum wavefront 
was 43; resequencing by one computer cycle reduced 
the initial maximum wavefront to 30. In another test, 
the initial sequencing was scrambled at random, 
which produced a maximum wavefront of 92; one 
computer cycle reduced this to 40. Subsequently, as 
the result of six additional computer cycles, the max- 
imum wavefront was reduced to 27. The computation 
time per cycle varies considerably, depending upon 
how quickly the connectivity matrix fills or upon the 
number of nodes that are processed before it can be 
determined that the current cycle is to be aborted 
because there will be no reduction of an existing 
maximum wavefront. 

Experience with the computer program indicates 
that minimum growth sequencing is effective, rapid, 
and capable of producing economies, such as in the 
reduction or elimination of lengthy effort to study a 
particular structural model for the purpose of gener- 
ating an acceptable sequencing pattern, and also in 
the reduction of computation time as the result of 
processing a relatively compact stiffness matrix. 
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